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detailed  simulation  of  the  radar  detection  algorithms  is  necessary  to  evaluate  a  radar  design  strategy  to 
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SECTION  1 


INTRODUCTION 


The  natural  ionosphere  can  produce  a  variety  of  disturbances  to  radar  sig¬ 
nal  propagation.  Disturbances  due  to  mean  or  very  large-scale  ionization  include 
attenuation,  phase  shift,  time  delay,  dispersion,  polarization  rotation,  refraction,  and 
multi-path.  In  addition,  relatively  small-scale  ionospheric  structure  in  the  propagation 
medium  can  cause  signal  scintillation-essentially  random  fluctuations  in  the  received 
signal  phase,  amplitude,  angle-of-arrival  and  other  signal  properties.  The  effects  of 
mean  propagation  disturbances  have  been  the  subject  of  many  studies  [ Budden ,  1985; 
Lawrence,  et  al .,  1964]  and  are  well  known.  The  effects  of  scintillation  on  the  detection 
performance  of  space  based  radar  are  the  subject  of  this  report. 

Electron  density  structure  in  the  natural  ionosphere  can  produce  random 
variations  in  the  amplitude  and  phase  of  a  propagating  wave,  even  at  frequencies  in 
the  gigahertz  range  [ Skinner ,  et  al.,  1971;  Taur,  1976].  These  rapid  variations  in  signal 
amplitude,  phase,  and  angle-of-arrival  are  called  scintillations  and  are  often  observed 
over  satellite  links  through  the  ambient  ionosphere  at  VHF  and  UHF  [Fremouw  et 
al.,  1978].  Radar  measurements  taken  with  the  ALTAIR  radar  in  the  Pacific  Test 
Range  are  known  to  give  severely  disturbed  scintillation  at  156  and  415  MHz  [Towle, 
1980].  Strong  scintillation  is  occasionally  observed  at  frequencies  as  high  as  C-band 
[ Franke ,  et  al.,  1984].  Since  even  small  fluctuations  in  the  received  signal  can  degrade 
performance,  the  effect  of  scintillation  must  be  considered  in  the  design  of  a  space 
based  radar  system  intended  to  operate  through  an  ionospheric  channel. 

The  severity  of  the  fluctuations  depends  upon  the  irregularity  of  the  ion¬ 
ization  structure  and  on  the  radar  geometry  and  frequency.  If  the  propagation  envi¬ 
ronment  is  highly  disturbed  and  the  radar  frequency  not  sufficiently  high,  worst-case 
scintillation  may  occur  wherein  the  received  signal  quadrature  components  are  uncor¬ 
related  Gaussian  variates,  after  one-way  propagation  over  the  severely  disturbed  path. 
Worst-case,  i.e.,  Rayleigh  amplitude,  scintillation  is  likely  to  occur  if  the  ionosphere 
is  highly  disturbed,  as  for  example  by  high  altitude  nuclear  explosions  [ Arendt  and 
Svicher,  1964;  King  and  Fleming,  1980]  or  by  chemical  releases  [  Wolcott,  et  al.,  1978]. 
An  increase  in  the  radar  frequency  or  a  lessening  of  the  ionization  irregularity  can  lead 
to  a  decreased  disturbance  in  the  received  signal  with  a  corresponding  change  in  the 
signal  statistical  description.  In  general,  however,  a  signal  at  one  frequency  (or  time) 
may  show  some  statistical  correlation  to  a  signal  received  at  a  different  frequency  (or 
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time).  The  effects  of  signal  decorrelation  with  time  and  frequency  at  all  levels  of  scin¬ 
tillation  severity  are  important  to  any  UHF  through  X-band  radar  system  that  must 
operate  through  an  ambient  or  disturbed  ionospheric  channel. 

Previous  authors  ( Schwartz ,  1956;  Linder  and  Swerling,  1956;  Walker,  1971] 
have  studied  the  performance  of  M  out  of  N  radar  detection  for  the  case  of  an  undis¬ 
turbed  propagation  channel.  In  an  earlier  work  [ Dana  and  Knepp,  1983],  the  impact 
of  strong  (Rayleigh)  fading  on  radar  noncoherent  detection  performance  was  reported, 
assuming  various  radar  geometries  and  target  models.  The  effects  of  frequency  corre¬ 
lation  were  also  studied  in  this  work,  however,  the  method  used  to  model  frequency 
correlation  was  limited  to  a  single  coherence  function. 

The  effects  of  scintillation  on  monostatic  radar  employing  “double-threshold” 
(Af  out  of  N)  detection  or  noncoherent  integration  are  the  subjects  of  this  work. 
(In  general,  noncoherent  integration  provides  superior  performance;  however,  double- 
threshold  detection  is  of  interest  since  it  affords  a  radar  design  simplification.)  In 
particular,  methods  are  developed  to  study  the  effects  of  varying  levels  of  scintillation 
on  the  probability  of  detecting  targets  of  constant  cross  section  as  well  as  targets 
whose  cross  section  follows  the  Swerling  II  model.  In  addition,  simulation  methods 
are  developed  to  study  the  effects  on  probability  of  detection  of  correlation  of  the 
propagation  channel  over  time  and  frequency. 

This  report  utilizes  the  Nakagami-m  probability  distribution  \N akagami, 
I960:  to  describe  the  statistics  of  the  received  signal  power  for  the  monostatic  radar  ge¬ 
ometry.  This  simple,  one-parameter  distribution  correctly  describes  both  the  weak  and 
strong  scattering  limits  of  the  actual  probability  distribution  of  scintillation  [ Knepp 
and  Valley,  1978]  and  has  been  found  to  give  a  very  good  fit  to  observed  scintilla¬ 
tion  data  I Fremouw ,  et  al.,  1980]  over  a  wide  range  of  scintillation  conditions.  The 
Nukagami-m  distribution,  as  used  here,  permits  the  study  of  radar  performance  over 
the  entire  range  of  possible  scintillation  conditions. 

When  the  two  extreme  assumptions  on  correlation  of  propagation  effects  are 
made  (i.e.,  receive  signal  statistics  due  to  propagation  fluctuations  for  different  bursts 
within  a  look  are  either  independent  or  else  completely  correlated),  then  expressions 
may  be  developed  for  detection  probability  which  are  readily  evaluated  via  numerical 
methods.  For  the  case  of  partially  correlated  bursts,  Monte  Carlo  methods  are  required 
to  determine  the  effect  on  probability  of  detection.  The  simulation  methods  developed 
here  allow  a  completely  general  specification  of  burst-to-burst  correlation  and  may  be 
used  to  study  the  effects  of  scintillation  on  the  performance  of  frequency  agile  radar. 
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In  this  report  a  coherent  pulse  train  is  referred  to  as  a  burst,  and  a  number  of 
bursts  that  are  combined  noncoherently  in  some  fashion  is  called  a  look.  It  is  assumed 
that  there  is  no  signal  decorrelation  over  the  duration  of  a  burst  and,  therefore,  no 
coherent  integration  loss  [ Dana  and  Knepp ,  1983]. 


1.1  IONIZATION  IRREGULARITY  DESCRIPTION. 


Figure  1  depicts  the  geographic  distribution  of  ionospheric  scintillation  as 
it  is  currently  known  (/larons,  1975],  The  severity  of  the  scintillation  is  indicated  in 
the  figure  by  the  density  of  the  cross-hatching.  Most  severe  ionization  irregularities 
occur  in  a  geographic  region  of  about  20  degrees  of  latitude  north  and  south  of  the 
geomagnetic  equator.  Equatorial  scintillation  is  observed  during  a  period  of  approxi¬ 
mately  8  PM  local  time  to  about  2  AM  and  may  regularly  take  the  form  of  saturated 
Rayleigh  amplitude  fading  at  VHF  and  UIIF,  and  occasionally  at  L-band.  Equatorial 
scintillation  is  known  to  be  more  severe  than  mid-  or  high-latitude  scintillation. 

At  high  latitudes,  scintillation  occurs  within  the  auroral  region  and  over 
the  polar  cap.  At  present  the  occurrence  of  auroral  scintillation  is  less  understood 
than  equatorial  scintillation.  However,  it  is  known  that  auroral  scintillation  is  more 
irregular  than  equatorial  scintillation  and  can  occur  anytime  during  the  local  day  or 
night.  Increased  magnetic  activity  brings  higher  levels  of  scintillation  and  can  cause 
changes  in  the  location  of  scintillation  boundaries  at  high  latitudes.  The  severity  of 
scintillation  can  be  measured  by  the  S4  scintillation  index  which  characterizes  the 
depth  of  fading. 

An  example  of  the  hourly  and  seasonal  distribution  of  amplitude  scintillation 
is  shown  in  Figure  2  [ Hawkins  and  Mullen,  1974]  in  terms  of  average  contours  of  the 
AFGL  SI  scintillation  index.  This  measurement  of  scintillation  severity  can  easily  be 
obtained  from  strip  charts.  The  innermost  contour  shown  corresponds  to  a  value  of  SI 
of  70  or  roughly  to  a  value  of  S4  of  0.36  signifying  8  dB  peak-to-peak  fading.  This  data 
was  obtained  from  an  accumulation  of  six  years  of  observations  of  ATS-3  from  Uuan- 
cayo,  Peru.  The  figure  shows  that  scintillation  at  this  location  is  most  severe  during 
the  months  of  February  through  October  and  during  the  hours  immediately  preceding 
local  midnight.  Similar  measurements  from  other  geographic  locations  indicate  the 
presence  of  a  longitudinal  dependence  on  seasonal  scintillation  activity.  Although  av¬ 
erage  results  of  this  type  are  extremely  useful  they  are  insufficient  to  design  an  SBll 
to  operate  during  signal  scintillation.  SBR  design  and  evaluation  are  best  served  with 
measurements  of  the  distribution  function  of  54  so  that  the  designer  can  determine 
the  percentage  of  time  peak-to-peak  signal  fluctuations  of  a  certain  level  are  exceeded. 
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depth  of  scintillation  FADING  (  proportional  to 
DENSITY  OF  CROSSHATCHING) 


Figure  1.  Geographic  distribution  of  ionospheric  fading. 


Fign  e  2.  Average  scintillation  (SI)  contours  at  136  MHz  obtained  over 
a  six  year  period.  Universal  time  shown  is  5  hours  less  than 
local  time  (05  UT  is  24  LT). 
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The  DNA  Wideband  satellite  experiment  has  assembled  the  most  extensive 
collection  of  phase  coherent  scintillation  data  to  date  | Fremouw  et  a/.,  1978].  Fig¬ 
ure  3  j Livingston ,  1978]  shows  the  seasonal  dependence  of  scintillation  at  137  MHz  as 
measured  during  Wideband  passes  observed  from  Ancon,  Peru.  The  Wideband  orbit 
is  sun-synchronous  so  that  the  satellite  always  passes  over  the  equator  at  a  constant 
local  time  of  about  11:30  PM.  Thus  from  a  fixed  ground  location,  three  or  four  passes 
are  observable  each  night  of  operation.  In  addition,  each  pass  takes  about  15  minutes. 
For  the  purposes  of  data  analysis  this  time  period  is  divided  into  many  10  second 
segments.  The  S4  scintillation  index  is  obtained  once  for  each  10  second  segment 
giving  many  measurements  of  S<  during  the  night.  Figure  3  shows  the  percentage  of 
measurements  of  at  Ancon  where  the  value  of  0.3  was  exceeded  at  137  MHz.  This 
value  corresponds  to  about  6  dB  peak-to-peak  fading. 


**'«J  -■  ?Ov  ZjZ  i  jAf.  f  £  &  N'AP\  Asnr.  ?.'A>  .'Uh: 

CALENDAR  MONTH 


Figure  3.  Seasonal  dependence  of  scintillation  severity  at  Ancon,  Peru. 

Values  of  scintillation  index  depend  on  geometry,  frequency,  and  the  ion¬ 
ization  irregularity  structure.  Results  from  the  DNA  Wideband  satellite  are  very 
appropriate  for  the  design  of  systems  that  might  use  high  inclination  orbits  since  the 
Wideband  satellite  was  in  an  almost  polar  orbit  at  an  altitude  of  1000  km.  Therefore 
these  measurements  may  be  used  directly  for  satellites  in  similar  orbits,  but  need  to 
be  scaled  to  apply  to  different  radar  transmission  frequencies  and  satellite  altitudes. 
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Figures  4  and  5  [ Knepp ,  1978]  show  Wideband  data  taken  on  a  particularly 
severe  pass  over  Ancon,  Peru  on  December  16,  1976.  Figure  4  shows  values  of  the  S4 
scintillation  index  obtained  at  138,  379  and  447  MHz  during  a  brief  portion  of  a  satel¬ 
lite  pass.  In  this  study  a  measurement  of  S4  is  obtained  by  averaging  over  each  10.5 
seconds  of  data  giving  sixteen  measurements  at  each  of  the  three  frequencies  during 
the  2.8  minute  time  period  shown.  During  the  beginning  of  this  period  propagation 
conditions  were  quite  severe  with  worst  case  Rayleigh  fading  (unity  S4)  observed  at 
VHF  and  the  two  UHF  frequencies  shown.  However  a  few  minutes  later,  the  scintilla¬ 
tion  severity  is  observed  to  decrease  somewhat  with  very  little  scintillation  (small  S4) 
at  UHF  at  23:51  although  at  VHF  scintillation  is  still  severe. 

Figure  5  shows  values  of  the  cross-correlation  function  of  the  received  inten¬ 
sity  or  power  as  obtained  from  the  first  six  10.5  second  intervals  shown  in  the  previous 
figure.  The  cross-correlation  coefficient  is  shown  for  correlations  between  the  lowest 
Wideband  UHF  tone  at  379  MHz  and  the  tones  at  390,  413,  436  and  447  MHz.  From 
the  figure  it  is  evident  that  the  tones  are  well  decorrelated  for  the  two  earliest  measure¬ 
ments  at  23:48:29  and  23:48:50  for  the  379-413  MHz  frequency  pair.  For  later  times 
during  this  Wideband  pass  the  UHF  frequency  tones  become  well  correlated  across 
the  spectrum  from  379  to  447  MHz  as  shown  by  the  curve  at  23:49:11.  The  data 
here  show  an  example  where  the  tones  at  379  and  413  MHz  were  decorrelated.  The 
coherence  bandwidth  is  proportional  to  transmission  frequency  to  the  fourth  power  so 
that  if  this  34  MHz  frequency  separation  were  scaled  to  apply  to  a  VHF  radar  at  200 
MHz,  the  resulting  coherence  bandwidth  would  be  about  4  MHz.  Thus  during  the 
scintillation  conditions  observed  during  these  measurements,  the  propagation  channel 
would  effectively  decorrelate  for  VHF  frequencies  separated  by  more  than  4  MHz.  A 
good  radar  design  could  take  advantage  of  this  decorrelation  with  transmission  fre¬ 
quency  by  noncoherently  combining  the  radar  returns  from  independent  samples  of 
the  fading  channel. 

The  mechanisms  accounting  for  irregularity  production  have  been  the  sub¬ 
ject  of  many  research  papers  over  the  last  decade  and  are  now  reasonably  well  un¬ 
derstood  ( Ktskintn  and  Ossakow ,  1981;  Tsanoda ,  1988] .  In  the  equatorial  regions  the 
dominant  mechanism  for  production  of  large  scale  structure  or  equatorial  spread-F 
is  the  collisional  Rayleigh-Taylor  instability.  The  emerging  view  is  that  high-latitude 
irregularities  in  electron  density  are  produced  primarily  by  the  gradient-drift  (or  E  x 
B)  instability. 

Although  other  forms  have  been  observed,  the  power  spectral  density  of 
the  electron-density  structure  is  well  represented  by  a  two-component  power-law  form 
with  the  irregularities  highly  elongated  along  the  direction  of  the  earth’s  magnetic 
field.  The  precise  values  of  the  parameters  that  describe  the  spectrum  are  under 
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Figure  4.  Simultaneous  VHF  and  UHF  scintillation  observed  during  the 
DNA  Wideband  satellite  experiment. 


Figure  5.  Cross-correlation  of  power  at  different  UHF  frequencies  during 
the  DNA  Wideband  satellite  experiment. 
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current  investigation  by  a  number  of  researchers.  However,  the  current  belief  is  that 
the  spectrum  behaves  as  K~*  for  spatial  wavenumbers  (K)  greater  than  the  freezing 
scale  and  flattens  to  between  K~2  s  and  K~2  7  for  scale  sizes  ranging  from  the  freezing 
scale  to  the  outer  scale.  The  outer  scale  ranges  from  10-100  km  and  the  freezing  scale 
ranges  from  150-1000  m.  One  view  of  the  irregularity  structure  causing  scintillation 
is  a  thick  layer  (~200  kilometers)  extending  upwards  from  the  base  of  the  F-region 
and  having  electron  density  fluctuations  on  the  order  of  50-100  percent  of  the  mean 
electron  density  [Basu  and  Basu,  1976], 
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SECTION  2 


f 

RECEIVED  SIGNAL  DESCRIPTION 


If  the  radar  signal  passes  through  a  disturbed  propagation  channel  subject 
to  electron  density  irregularities,  the  received  signal  amplitude  may  experience  fluc¬ 
tuations  or  scintillation.  Knowledge  of  the  probability  distribution  of  the  received 
signal  is  essential  to  determine  radar  system  performance.  Since  the  1950s,  many 
authors  have  investigated  the  probability  distribution  of  the  received  signal  after  one¬ 
way  propagation  through  turbulence  [ Knepp  and  Valley,  1978  and  references  therein]. 
A  large  number  of  probability  distribution  functions  that  describe  various  aspects  of 
the  received  scintillated  signal  are  possible.  Many  distributions  correctly  describe  the 
weak-scattering  regime,  and  parametrically  extend  to  the  strong-scattering  limit.  The 
joint  Gaussian  distribution  is  generally  the  most  useful  to  describe  higher  moments 
of  the  received  electromagnetic  field.  However,  in  a  definitive  paper  Fremouw,  et 
al.  [1980]  utilized  measurements  to  show  that  the  Nakagami-m  distribution  is  clearly 
superior  to  the  log-normal,  the  generalized  Gaussian  and  the  two-component  Gaussian 
distribution  to  describe  the  probability  density  function  of  the  received  power  for  the 
case  of  radio  wave  scattering.  This  conclusion  is  based  on  a  comprehensive  analysis  of 
the  large  body  of  multi-spectral  scintillation  data  collected  during  the  DNA  Wideband 
Satellite  experiment. 


2.1  RECEIVED  SIGNAL  FIRST-ORDER  STATISTICS. 


In  fading  conditions  it  is  convenient  to  write  the  received  signal  power  as 
[ Dana  and  Knepp,  1983] 


ST  —  S0Sa/(a)  (1) 

In  this  expression  for  the  received  power,  50  is  the  mean  signal  power  received  from  the 
target,  S  is  the  fractional  change  in  the  signal  power  due  to  variations  caused  by  the 
propagation  channel,  and  o/(a )  is  the  fractional  change  in  the  signal  power  caused 
by  target  cross  section  fluctuations.  The  factor  So  contains  the  mean  signal  level; 
therefore,  the  mean  value  (S)  is  unity  Since  the  mean  signal  level  has  no  statistical 
variation,  the  received  signal  power  fluctuations  can  be  expressed  as  the  product  of 
the  effects  of  fluctuations  due  to  propagation  disturbances  and  to  target  variations. 
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2.2  TARGET  STATISTICS. 


Consider  the  case  of  a  radar  that  transmits  bursts  that  are  separated  by 
a  iarge  amount  of  time  relative  to  the  decorrelation  time  of  the  target  cross  section 
variations.  A  similar  situation  arises  for  the  case  that  the  radar  transmission  fre¬ 
quency  is  changed  sufficiently  from  burst  to  burst,  so  that  the  target  cross  section  is 
decorrelated  from  burst  to  burst.  In  either  case  the  target  cross  section  variation  can 
often  be  assumed  to  be  described  by  a  Swerling  II  model.  A  Swerling  II  model  ap¬ 
plies  to  the  case  that  the  cross  section  varies  independently  from  burst  to  burst.  This 
nomenclature  differs  from  the  original  Swerling  II  convention  only  in  the  replacement 
of  Swerling’s  “pulse”  by  a  “burst”  in  this  work.  The  probability  density  function  for 
target  cross  section  variation  is  given  by  the  expression 

p{°)  =  -Ayexp(-a/(o»,  a  >0  (2) 

where  (o)  is  the  mean  value.  The  angle  brackets  denote  a  stochastic  average.  Since  a 
target  with  constant  cross  section  is  also  of  interest  to  certai  r  defense  radars,  results 
are  also  included  for  this  case. 


2.3  PROPAGATION  CHANNEL  STATISTICS. 


To  completely  characterize  the  first-order  statistics  at  the  radar  receiver, 
the  statistics  of  the  signal  fluctuations  caused  by  propagation  through  a  disturbed 
channel  are  required.  Using  Nakagami-m  statistics  the  probability  density  function  of 
fading  on  the  one-way  propagation  channel  is  controlled  by  the  m-parameter  defined 
by  the  equation 


2  i(P-(P))2) 

4  (py 


(3) 


In  this  equation  SA  is  referred  to  as  the  scintillation  index  and  is  a  measure  of  the 
severity  of  the  fluctuations  of  power.  The  quantity  P  is  the  received  power  on  a  one¬ 
way  propagation  path.  The  S <  index  is  thus  the  normalized  standard  deviation  of 
the  received  power  on  a  one-way  propagation  path  where  the  transmitted  power  is 
constant.  Values  of  S+  generally  range  from  a  minimum  of  zero  signifying  constant 
power  or  no  scintillation  to  a  maximum  of  unity  indicating  worst-case  (Rayleigh) 
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fading  of  amplitude  where  the  in-phase  and  quadrature  components  of  the  received 
one-way  signal  are  uncorrelated  Gaussian  variates.  Values  of  S+  greater  than  unity 
have  been  observed,  but  these  indicate  the  presence  of  signal  enhancements  due  to 
focusing.  Since  S4  has  been  used  more  frequently  in  the  past  to  describe  scintillation, 
it  will  be  used  to  quantify  the  scintillation  severity  in  place  of  m  in  the  remainder  of 
this  report. 

For  a  one-way  propagation  path  the  probability  density  function  of  received 
power  is  given  by  [Nakagami,  1960] 


mmSm~1  f  —mS )  . 

Pl (s)  =  r(m)(s)» exp  V ~w }  ’  -  ’  one~way  (4) 

where  the  subscript  refers  to  the  one-way  propagation  path.  For  the  case  of  monostatic 
radar  operation,  the  transmitter  and  receiver  are  collocated  so  the  signal  propagates 
twice  over  the  same  path,  passing  through  identical  irregularities.  In  this  case  the  re¬ 
ceived  voltage  is  proportional  to  the  square  of  the  voltage  after  one-way  propagation 
[ Knepp ,  1985].  The  probability  density  function  for  the  received  power  may  be  ob¬ 
tained  from  Equation  4  using  the  transformation  Q(monostatic)  =  S'2 (one-way)  with 
the  result 


[m(m  -f  i)\ml2Qml2  1 
MQ)  =  jr(m)  W)W.  exp  < 


m{m  +  1  )Q 

N  W) 


,Q>o,  two-way  (5) 


Note  that  the  Nakagami-m  density  is  a  generalization  of  the  chi-square  prob¬ 
ability  density  function.  The  chi-square  density  can  be  obtained  from  Equation  4  by 
replacing  the  m-parameter  with  nf  2  and  performing  the  change  of  variable  y  =  nS 
on  the  resulting  equation.  In  the  resulting  equation,  y  is  the  chi-square  variate  and  n 
is  the  number  of  degrees  of  freedom  which  is  restricted  to  positive  integer  values. 

It  is  straightforward  to  compute  the  relationship  between  the  value  of  S4 
measured  on  the  monostatic  radar,  two-way  path  as  a  function  of  m  on  the  one-way 
path  as  ^(two-way)  =  (4m2  +  10m  +  6)/(m(m  + 1)2).  Thus,  under  Rayleigh  statistics, 
the  maximum  value  of  S4  on  the  two-way  path,  corresponding  to  a  monostatic  radar 
geometry,  is  y/E.  Figure  6  shows  the  cumulative  distribution  of  the  received  power  on  a 
round  trip  propagation  path  that  characterizes  a  monostatic  radar  geometry  consisting 
of  two  one-way  propagation  paths.  The  cumulative  distribution  is  shown  as  a  function 
of  the  value  of  S4  on  each  one-way  propagation  path.  For  worst-case  scintillation  in 
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the  monostatic  radar  geometry  it  is  seen  from  the  figure  that  the  probability  of  a 
10  dB  fade  or  greater  is  about  40  percent.  This  high  probability  of  fading  is  reflected 
in  degradation  of  various  aspects  of  radar  performance  including  reduced  radar  target 
detection  performance. 


oi  i  '  c  *c  e.c  ss  i  a  i~- 


PROBABILITY  [POWER  \FCWER S  .V  ] 

Figure  6.  Cumulative  probability  distribution  of  received  two-way  power 
as  a  function  of  on  the  one-way  propagation  path. 


2.4  PROPAGATION  CHANNEL  COHERENCE. 


The  above  discussion  fully  describes  the  first-order  power  statistics  to  be 
expected  after  propagation  of  a  radar  signal  through  a  disturbed  ionospheric  channel. 

The  second-order  fading  statistics  are  specified  by  the  correlation  function 
of  the  received  complex  voltage.  For  the  case  of  one-way  propagation  of  an  initially 
constant  amplitude  signal  through  a  severely  disturbed  ionospheric  channel,  the  auto¬ 
correlation  function  of  the  received  voltage  is  given  a s  the  two-position,  two-frequency 
mutual  coherence  function.  The  effective  velocity  of  the  line  of  sight  of  the  radar  signal 
through  the  ionospheric  irregularities  can  be  utilized  to  convert  the  spatial  coordinates 
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of  the  mutual  coherence  function  into  temporal  coordinates  and  thereby  obtain  the 
correlation  function  of  signal  fluctuations  due  to  ionospheric  irregularities. 

For  worst-case  Rayleigh  fading,  the  correlation  function  of  the  received  com¬ 
plex  voltage  E(t,f)  has  the  form 

(B{t  +  r,f  +  /,)£•(<,/)>  =  l/Tot3  exp  {_ J  (,  +  ,/,//„,)-■  (6) 

where  r0  is  the  decorrelation  time  (the  fading  rate  is  1/tq)  for  fluctuations  over  a  one¬ 
way  propagation  path.  The  actual  value  of  r0  is  a  function  of  radar  geometry  and  of 
the  irregularity  structure  and  intensity  of  the  disturbed  ionospheric  channel.  Large 
values  of  tu  correspond  to  slow  fading  conditions  and  small  values  correspond  to  fast 
fading. 

As  a  concrete  example,  consider  the  case  of  a  radar  and  target  separated 
by  a  layer  of  ionization  as  might  occur  in  the  case  of  an  SBR  observing  a  target 
near  the  ground.  For  a  K~<  in  situ  power  spectrum  of  three-dimensional  ionization 
irregularities  between  outer  scale  L0  and  inner  scale  the  decorrelation  time  is 

r0  =  y/2L0/\Jln{L0/l,)o^vL  (7) 

where  vi  is  the  velocity  of  the  line  of  sight  through  the  center  of  the  ionized  layer,  o ^  is 
2(re A)2L0LA N?  rad2,  A  is  the  RF  wavelength,  rt  is  the  classical  electron  radius  (2.82 
x  10“ 15  m),  L  is  the  thickness  of  ionized  layer,  and  AN?  is  the  variance  of  electron 
density  irregularities. 

It  is  assumed  that  r0  is  large  with  respect  to  the  duration  of  the  transmitted 
pulse.  The  received  signal  is  then  coherent  during  the  pulse  duration  that  is  typically 
of  the  order  of  several  tens  of  microseconds. 

The  channel  coherence  bandwidth  is  a  measure  of  the  maximum  bandwidth 
available  in  the  propagation  channel  over  which  it  is  possible  to  transmit  a  signal 
without  undesired  pulse  distortion.  That  is,  in  a  fading  environment,  signal  spectral 
components  separated  by  less  than  the  coherence  bandwidth  exhibit  correlated  fluc¬ 
tuations.  If  the  signal  spectral  components  are  separated  in  frequency  by  an  amount 
greater  than  the  coherence  bandwidth,  different  spectral  components  will  undergo 
uncorrelated  fading.  This  distortion  in  the  received  signal  spectrum  causes  the  re¬ 
ceived  time  domain  signal  to  display  undesired  time  sidelobes.  On  the  other  hand,  a 
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propagating  pulse  remains  undistorted  as  long  as  the  maximum  instantaneous  signal 
bandwidth  is  less  than  the  coherence  bandwidth.  This  is  often  the  case  for  frequency- 
agile  radar  systems.  In  addition,  frequency  agile  radar  signals  that  are  separated  in 
frequency  by  an  amount  exceeding  the  coherence  bandwidth  will  experience  nearly 
independent  fading.  In  this  work  it  is  assumed  that  the  channel  coherence  bandwidth 
is  sufficiently  large  compared  with  the  radar  pulse  bandwidth  so  that  no  time-domain 
distortion  of  the  received  pulse  occurs. 

For  a  in  situ  power  spectrum  of  three-dimensional  ionization  irregu¬ 

larities  and  one-way  signal  propagation  path  geometry  with  the  radar  and  target  on 
opposite  sides  of  a  scattering  layer,  the  coherence  bandwidth  is  given  by 


,  _ tcc(z,  zr)L0 _ 

~  T]\*\n(L0;it)ztzrL'KN} 

In  this  expression,  c  is  the  velocity  of  light  in  vacuum,  zt  is  the  distance  from  the 
transmitter  to  the  center  of  the  ionized  layer,  zr  is  the  distance  from  the  target  to 
the  center  of  the  ionized  layer,  and  zt  +  zr  is  the  total  one-way  propagation  distance. 
The  relationship  between  one-  and  two-way  propagation  paths  is  discussed  in  Knepp , 
[  1985 1 .  To  obtain  the  appropriate  description  of  a  two-way  monostatic  radar  propa¬ 
gation  geometry,  reduce  the  values  of  r0  and  fcoh  on  the  one-way  path  by  a  factor  of 
v/2. 


(S) 


In  general  the  signal  power  fluctuation  due  to  variations  caused  by  the  prop¬ 
agation  channel,  5,  exhibits  correlation  over  both  time  and  frequency  as  discussed 
above.  However,  in  this  work  only  the  burst-to-burst  correlation  is  of  interest.  There¬ 
fore,  the  resulting  correlation  function  may  be  modeled  as  a  frequency  coherence 
function  given  by 


Css(fd)  —  (S{f)S(f  +  /<*))  -  (S)2  (0) 

where  it  is  assumed  that  S  is  a  stationary  process  over  the  frequency  range  of  interest. 

Consideration  of  channel  coherence  is  important  for  the  assessment  of  the 
effectiveness  of  burst  combining  using  a  frequency  agile  radar  in  a  fading  channel. 
If  the  separation  of  burst  transmission  frequencies  is  great  enough  to  insure  that 
S(f)  is  uncorrelated,  then  radar  performance  will  be  greatly  enhanced  through  burst 
combining.  However,  if  the  received  power  shows  correlation  from  burst  to  burst,  then 
it  is  important  to  account  for  Css  in  the  calculation  of  detection  probability. 
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SECTION  3 


RADAR  SYSTEM  CHARACTERISTICS 


In  this  report  it  is  assumed  that  an  SBR  is  required  to  operate  through  a 
disturbed  ionospheric  channel  to  detect  and  track  targets  near  the  earth’s  surface. 
Thus,  relative  to  a  ground  based  radar  with  similar  functions,  an  SBR  has  several 
limitations.  First,  targets  are  detected  and  tracked  at  very  long  ranges.  Second, 
available  onboard  transmitter  power  is  relatively  low.  Third,  because  of  the  large 
target  ranges  involved,  large  areas  of  the  earth’s  surface  are  illuminated  with  resultant 
large  clutter  returns,  even  with  a  very  narrow  antenna  beamwidth. 

The  first  two  points  imply  low  received  signal-to-noise  ratio  per  pulse  and 
therefore  require  long  integration  times.  However,  the  cross  section  of  a  moving  target 
remains  constant  or  coherent  for  only  a  few  tens  of  milliseconds  because  of  target  mo¬ 
tion  and  resulting  constructive  and  destructive  interference  between  many  scattering 
centers.  Hence,  during  a  radar  look,  the  total  energy  transmitted  at  a  target  is  divided 
into  a  number  of  bursts.  Each  burst  consists  of  some  number,  n,  of  pulses  which  are 
coherently  integrated.  This  group  of  pulses  will  be  referred  to  as  a  burst.  Multiple 
bursts  may  be  transmitted  to  form  a  look,  with  each  burst  transmitted  at  a  different 
frequency  in  the  case  of  a  frequency  hopping  radar.  The  detected  amplitude  of  all  the 
bursts  which  form  the  total  radar  target  look  are  then  noncoherently  combined  in  a 
postdetection  integration  process.  The  resulting  signal  power  is  then  compared  to  a 
threshold  level  to  decide  whether  or  not  a  target  has  been  detected  during  the  radar 
look. 


One  reason  for  choosing  a  waveform  consisting  of  many  coherently  related 
pulses  is  that  the  radar  transmitter  may  be  power  limited  and  unable  to  generate 
sufficient  energy  in  a  single  pulse.  Thus  many  pulses  with  low  power  are  coherently 
transmitted,  then  coherently  integrated  on  reception  to  achieve  high  signal-to-noise 
ratio.  A  second  advantage  in  using  a  coherent  pulsed  radar  waveform  is  that  doppler 
processing  techniques  may  be  applied  to  the  pulses  comprising  a  single  burst  to  achieve 
clutter  suppression. 

The  use  of  frequency  hopping  can  produce  independent  samples  of  the  tar¬ 
get  cross  section  from  burst  to  burst.  Depending  on  the  target  geometry  and  motion 
relative  to  the  radar,  the  target  cross  section  from  burst  to  burst  may  also  be  indepen¬ 
dent.  In  either  case  a  Swerling  II  cross  section  model  applies.  Since  a  nonfluctuating 
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target  serves  as  a  convenient  reference  for  computation  and  is  also  of  interest  for  many 
types  of  defense  systems,  results  are  also  given  for  this  target  model. 

When  the  time  duration  of  a  burst  is  short  enough  so  that  the  received 
signal  power  remains  essentially  constant  within  a  burst,  then  the  fading  is  referred 
to  as  slow  fading.  The  slow  fading  assumption  will  be  assumed  to  apply  for  all  cases 
considered  in  this  report. 

Two  types  of  radar  detection  techniques  are  considered  here,  namely  double 
threshold  detection  and  noncoherent  integration.  In  both  techniques  the  first  stage  is 
quadrature  detection  of  the  radar  return  from  a  transmitted  burst.  The  quadrature 
detector  output  is  a  measure  of  received  power.  For  the  M  out  of  N  detection  scheme, 
N  bursts  are  transmitted  to  form  a  look.  If  the  power  received  from  M  or  more 
bursts  exceeds  a  predetermined  threshold,  then  a  detection  is  declared.  In  the  case 
of  noncoherent  integration,  the  power  measured  from  each  burst  in  a  look  is  summed 
and  if  the  resultant  sum  exceeds  a  threshold,  then  a  detection  is  declared. 

A  phased-array  radar  system  having  sufficient  beam  pointing  flexibility 
might  utilize  a  modified  M  out  of  N  detection  scheme.  For  example,  if  a  it  is  declared 
on  the  first  M  bursts,  the  rest  need  not  be  transmitted.  For  a  radar  that  must  detec* 
multiple  targets,  additional  bursts  may  be  transmitted  to  a  specific  location  upon  a 
threshold  crossing  of  the  return  from  the  first  burst. 
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SECTION  4 


PROBABILITY  OF  DETECTION 


The  goal  of  this  section  is  to  develop  methods  to  obtain  the  probability 
of  detection  for  a  radar  operating  in  the  presence  of  scintillation.  Mathematical  ex¬ 
pressions,  which  are  readily  evaluated  using  numerical  methods,  are  found  for  the 
detection  probability  when  constraints  are  placed  on  the  burst  correlation  function 
of  the  random  process  5,  -  the  fractional  change  in  signal  power  in  burst  t  due  to 
variations  caused  by  the  propagation  channel.  For  the  remaining  cases,  the  detection 
probability  is  calculated  via  Monie  Carlo  simulation  methods.  The  simulation  method 
presented  here  is  general  enough  to  include  any  coherence  function  for  S. 

Expressions  that  may  be  readily  evaluated  numerically  are  developed  for  the 
following  cases:  M  out  of  N  detection  of  both  constant  and  Swerling  II  targets  when 
the  received  power  S  is  assumed  independent  from  burst  to  burst,  M  out  of  N  for 
both  constant  and  Swerling  II  targets  when  5  is  assumed  identical  for  each  burst  in  a 
look,  and  noncoherent  integration  of  bursts  for  a  Swerling  II  target  when  S  is  assumed 
identical  for  each  burst  in  a  look.  In  all  cases  5  is  assumed  to  follow  a  Nakagami-m 
amplitude  distribution. 

Monte  Carlo  methods  developed  here  may  be  used  to  determine  the  prob¬ 
ability  of  detection  for  both  M  out  of  N  detection  and  noncoherent  integration  for 
a  Swerling  II  target.  The  received  signal  power  fluctuation  due  to  propagation  for 
each  burst,  Si,  follows  an  arbitrary  correlation  function  and  a  Nakagami-m  amplitude 
distribution. 


4.1  RECEIVER  MODEL. 


Under  slow  fading  conditions  let  the  signal  amplitude  and  phase  during  a 
received  radar  burst  be  constant  so  that  the  in-phase  and  quadrature  voltage  compo¬ 
nents  are  given  by  the  expressions 


i  —  a  cos  <f>  +  n,  (10) 

q  =  asin  <f>  I  nq  (ll) 
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where  »  is  the  voltage  (signal  plus  noise)  of  a  single  burst  in  the  in-phase  channel  and 
q  is  the  voltage  of  a  single  burst  in  the  quadrature  channel.  The  total  average  signal 
power  is  (i2)  +  ( q 2)  —  Sr.  The  noise  is  assumed  to  be  additive  white  Gaussian  noise 
with  probability  density  function 


P(n.)  = 


\/2ttoh 


exp 


\  20  N 


(12) 


with  a  similar  expression  for  p{nq).  The  total  noise  power  per  burst  is 


K2)  +  («J)  =  2<r*  (13) 

so  that  the  signal-to-noise  ratio  per  burst  is  Sr/2a^.  Now  the  output  voltage  from 
a  quadrature  detector  giving  w  —  (»2  +  q2)1/2  has  the  well-known  Rician  probability 
density  function  conditioned  on  Sr 


PMSr)  =  —  exp  j-  ^  |  J0  (w^/Fr/aj^  (14) 

where  I0  is  the  modified  Bessel  function.  The  above  density  function  may  also  be 
conditioned  jointly  on  5  and  a  through  the  use  of  Equation  1.  For  a  Swcrling  II 
target,  the  probability  density  function  of  the  target  cross  section  is  given  by  Equation 
2.  Now  since  S  and  a  are  independent  variables,  a  conditional  density  with  respect 
to  5  alone  is  obtained  via 


roo 

p{w\S)=  p(o)  p{w\Sr  ~  \S0Sc/{o)})  do  (15) 

Jo 

where  Sr  =  SqSo/(o).  The  integral  expression  in  Equation  15  easily  reduces  to  the 
conditional  density  for  vu. 


p{w\S)  = 


w 


(1  +  SoS/20y 


exp 


-ur 


2  Oy  (1 


S0S/2a%) 


Swcrling  II  target  (10) 
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4.2  NONFLUCTUATING  TARGET. 


For  a  nonfluctuating  target  the  probability  density  function  of  the  target 
cross  section  can  be  represented  by  the  Dirac  delta  function  6  [o  —  S0).  In  this  case 
the  integration  in  Equation  15  gives 


p(to|5) --  ^-expj — |  /o  (w\/s0S /ax')  ,  constant  target  (17) 
°n  (.  *°n  j  '  ' 

Note  that  the  mean  received  power  for  both  the  Swerling  II  and  the  constant  target 
is  S05  -)-  2o2n  as  it  should  be. 


4.3  SINGLE  BURST  DETECTION. 


The  quadrature  detector  output  is  compared  to  a  preset  threshold  at  the 
end  of  a  hurst  interval,  and  a  “hit”  or  a  “miss”  is  declared.  The  threshold  is  set  on  the 
basis  of  the  noise  alone  so  that  the  probability  of  false  alarm  is  small  The  probability 
of  a  false  alarm  for  detection  of  a  single  burst  is  given  by 

P/a  =  /  p(u>|S  =  °)  dw  (18) 

Jth 

The  quantity  p]a  is  simply  the  single-burst  probability  that  the  noise  amplitude  alone 
exceeds  the  threshold.  The  above  equation  is  integrated  to  obtain  the  familiar  result 


th  -  y/-2a%  In  (p/o)  (19) 

Thus  the  threshold  amplitude  is  expressed  as  a  function  of  the  RMS  noise  amplitude 
and  the  desired  probability  of  false  alarm  for  a  single  burst.  The  threshold  is  the  same 
for  all  channel  and  target  models,  since  the  noise  is  independent  of  target  and  channel 
statistics. 


Equations  15  and  19  may  be  used  to  find  the  single-burst  probability  of 
detection  as  a  function  of  S ,  the  fractional  change  in  signal  amplitude  due  to  scintil¬ 
lation. 
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Pd[S)  -  f  p(w\S)dw 

Jth 

—  p(fl+sl'SNR'>'i  >  Swerling  II  target  or 

=  Q  ^-2\n{pfa),yj2S(SNR)')  ,  constant  target  (20) 


where  (SNR)  —  S0/2o2N  is  the  average  signal-to- noise  ratio  at  the  amplitude  detector 
output  and  Q  is  the  Marcum-Q  function  [Marcum,  1948].  When  5  =  1,  Equation  20 
is  the  familiar  formula  for  detection  probability  of  a  Swerling  II  target. 


4.4  M  OUT  OF  N  DETECTION. 


In  the  case  of  M  out  of  N  or  double  threshold  detection,  it  is  straightforward 
to  obtain  simple  analytic  expressions  for  the  overall  probability  of  detection  in  terms 
of  the  probability  of  detection  per  burst,  provided  that  the  radar  returns  from  each 
burst  are  independent. 

Now  assume  that  the  probability  of  detection  per  burst  is  given  by  pd .  To 
find  the  overall  probability  of  detection,  Pp,  le*.  Pb{*,j)  he  the  probability  of  obtain¬ 
ing  exactly  t  detections  in  j  bursts.  Then  the  probability  of  obtaining  M  or  more 
detections  out  of  N  bursts  is  the  sum 


p(M,N)  =  pb(M,N)  +Pb(M  +  l,N)  +  --+Pb(N,N) 

=  £  ft(*,AT)  (21) 

k=M 

but  since  pt,  is  given  by  the  binomial  distribution  [ Feller ,  1957], 

n(-\j)  =  (’)pi(i  -*)*-•'  (22) 

then  it  is  easy  to  obtain  the  probability  of  detection  for  the  M  out  of  N  process, 

Pv=  E  (!)  (p.)‘(l  -  P.r-*1  (23) 
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In  the  absence  of  a  target,  the  overall  probability  of  false  alarm  for  double 
threshold  detection  is  given  by  a  similar  expression 

=  £  (?)  (P/.)‘(1  -  Pt.r -*»  (24) 

k-M 

where  p/0  is  the  single-burst  false  alarm  probability. 

4.5  INDEPENDENT  BURSTS. 

Now  consider  the  case  that  the  contribution  to  the  received  signal  due  to 
channel  fluctuations  is  independent  from  burst  to  burst,  either  due  to  temporal  or 
frequency  offsets  between  bursts.  Here  the  probability  of  detection  for  a  single  burst 
is  obtained  by  averaging  pj(5),  given  by  Equation  20,  over  the  distribution  of  channel 
fluctuations  p2{S). 

P4=  f°°  P4(S)p2(S)dS  (25) 

JO 

where  p2{S)  is  given  by  Equation  5. 

4.6  CONSTANT  PROPAGATION  CHANNEL. 

Now  consider  the  case  that  the  contribution  to  the  received  signal  due  to 
the  disturbed  channel  is  constant  over  all  the  bursts  that  comprise  a  look.  In  this  case 
the  probability  of  detection  during  fading  is  computed  by  averaging  the  probability 
of  detection  for  constant  propagation  conditions  over  the  appropriate  distribution 
function  describing  the  propagation  channel  fluctuations.  In  this  case  the  overall 
detection  probability  may  be  written  as 

roo  f  N  . 

r»=  £  (?)  (MS))‘(1  -  P,(s)ds  (26) 

lk-M 

where  pj(S)  is  given  by  Equation  20. 
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4.7  NONCOHERENT  INTEGRATION  -  CONSTANT 
PROPAGATION  CHANNEL. 


For  a  Swerling  II  target  and  a  constant  propagation  channel,  Dana  and 
Knepp  [1986]  have  developed  a  general  formulation  for  the  case  of  noncoherent  in¬ 
tegration  that  may  be  applied  to  the  case  of  Nakagami-m  statistics.  The  detection 
probability  is  evaluated  using  the  expression 

Pd  =  I"  PD(S)p2(S)dS  (27) 

Jo 

where  p2(S)  is  given  by  Equation  5  and  where 


„  r{M,thHoUl  +  S{SNR))\) 

Pd(s]  =  r  (M) 

(28) 

r(a,I)  = 

(29) 

and  the  threshold  is  set  using 

PrA  =  T(M,th/<j2N)/T(M)  (30) 

r(m)  is  the  usual  gamma  function  of  a  single  argument  and  T(a,t)  is  an  incomplete 
gamma  function.  Equation  29  is  easily  inverted  numerically  to  give  thf  o7s  as  a  function 
of  the  probability  of  false  alarm  and  the  number  of  bursts  M. 


4.8  NUMERICAL  METHODS. 


It  is  straightforward  to  numerically  evaluate  the  foregoing  integral  expres¬ 
sions  for  probability  of  detection.  For  the  case  of  M  out  of  N  detection,  using  specified 
values  of  M,  N  and  the  overall  probability  of  false  alarm,  Pfa,  the  first  step  is  the  nu¬ 
merical  inversion  of  Equation  24  to  obtain  the  single-burst  probability  of  false  alarm, 
P/a.  If  desired,  the  threshold  is  then  obtained  using  Equation  19.  The  quantity  pja 
is  used  in  Equation  20  to  obtain  Pi[S)  and  the  integrals  of  Equations  25  and  26  are 
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performed  as  desired.  This  entire  process  is  repeated  as  often  as  desired  for  different 
values  of  signal-to-noise  ratio  and  scintillation  index. 

A  similar  procedure  is  used  for  the  case  of  the  noncoherent  combining  of 
bursts  from  a  Swerling  II  target  in  a  constant  propagation  channel. 


4.9  GENERAL  CORRELATION  OF  BURST  RETURNS. 


In  the  general  case  of  a  disturbed  propagation  channel,  the  radar  returns 
from  the  bursts  that  comprise  a  look  will  exhibit  some  degree  of  correlation  in  time  and 
frequency  as  discussed  earlier.  In  this  case,  the  use  of  numerical  integration  techniques 
becomes  prohibitive  because  the  probability  of  detection  is  expressed  as  multivariate 
integrals.  Therefore,  Monte  Carlo  simulation  techniques  must  be  used  to  obtain  the 
probability  of  detection  for  the  case  that  burst-to-burst  correlation  is  described  by 
a  general  correlation  function.  The  techniques  described  below  are  valid  for  general 
correlation  functions  of  any  form.  However,  in  the  examples  to  be  presented,  a  simple 
Gaussian  function  is  used.  This  simple  one-parameter  distribution  is  quite  useful  since 
a  Gaussian  function  is  often  used  to  fit  more  complicated  distributions. 

The  simulation  is  carried  out  by  generating  random  samples  of  the  received 
signal  power,  Sr,  for  each  burst  within  a  look.  The  resulting  samples  may  be  processed 
directly  using  either  double-threshold  detection  or  noncoherent  integration.  After  a 
sufficient  number  of  looks  are  processed,  an  estimate  is  obtained  of  the  detection 
probability  at  a  given  SNR  and  probability  of  false  alarm. 

In  order  to  have  a  valid  representation  of  the  fading  channel,  it  is  required 
that  the  statistics  of  received  power  follow  the  Nakagami-m  distribution.  Fortunately, 
as  noted  previously,  the  \2  number  of  degrees  of  freedom  is  related  to  the  Nakagami-m 
parameter  by  the  expression  m  =  n/2,  n  —  1,2,3  ■  • -.  A  useful  range  of  values  of  S4 
may  be  obtained  using  the  available  half-integer  values  of  m. 

A  variate  with  n  degrees  of  freedom  may  easily  be  generated  from  the 
expression 


S  =  v\  +  v\  +  •••  +  v2n  (31) 

where  the  v,  are  independent,  zero  mean  Gaussian  variates  with  variance  a l  selected 
so  that  the  mean  value  of  S  is  unity. 
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It  is  also  required  that  all  values  of  power  S  for  bursts  in  a  single  look  be 
related  according  to  the  arbitrary  correlation  function  Css ■  Write  5  as  a  function  of 
the  index  of  the  burst,  »,  where  there  are  N  bursts  per  look  and  t  —  0, 1,  •  •  • ,  N  —  1. 
Therefore,  it  is  required  that 

Cssti)  =  <S(»)S((»  +  i),  mod  N))  -  ( 5 )2 

J-  —  (TV  —  1),  •  •  • ,  (TV  —  1)  (32) 

where  the  notation  “mod”  refers  to  tiie  modulo  remainder  function  defined  by  (a,  mod  N) 
r  where  a  =  r  +  nN  and  n  is  any  integer  or  zero. 

Expanding  Equation  31  for  all  bursts  in  a  look  gives  a  set  of  equations 


5(1) 

=  v\{ 1) 

+ 

«!(i) 

+ 

®2(i) 

+  •  ■ 

■  •  + 

«*(!) 

5(2) 

-  "?(2) 

+ 

«I(2) 

+ 

vj(2) 

»*(2) 

(33) 

S(N) 

-  v\{N) 

+ 

vl(N) 

+ 

t'I(TV) 

+  •  ■ 

•  •  + 

Equation  33  may  be  satisfied  if 

Css  (j)  =  2nClu{j)  (34) 

where 

Cm  (j)  =  (vjk(i')u*  ((*  +  j),  mod  N));  any  k,  any  i  (35) 

^he  above  assertion  is  easily  proven  using  the  relation  [ Papoulis ,  1984] 

(*V)  =  (xJ)(y2)  +  2{xy)7  (36) 

where  x  and  y  are  jointly  normal  with  zero  mean. 

It  is,  therefore,  necessary  to  generate  n  gaussian  vectors,  each  of  length  TV, 
that  satisfy  Equation  34  that  are  to  be  used  in  Equation  33.  An  efficient  and  flexible 
method  for  generating  a  vector  of  jointly  normal  variates  is  presented  in  Hurst  and 
Knap  [  1972] . 
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SECTION  5 


RESULTS 


In  this  section  quantitative  results  are  given  for  probability  of  detection 
for  a  wide  range  of  scintillation  conditions  ranging  from  weak  to  strong  scattering 
and  for  various  amounts  of  decorrelation  between  bursts  in  a  look.  It  is  shown  that 
both  double  threshold  detection  and  noncoherent  combining  of  bursts  offer  a  degree  of 
mitigation  against  the  effects  of  scintillation  when  there  is  some  diversity  gain  available 
through  decorrelation  of  bursts  in  a  look.  Although  the  probability  of  false  alarm  is 
arbitrary  in  our  formulation  and  is  an  input  variable  to  the  simulations,  all  results  to 
be  presented  assume  that  the  detection  threshold  is  set  so  that  the  false  alarm  rate 
per  look  is  10~6. 

Figures  7  and  8  illustrate  the  effects  of  various  levels  of  fading  on  a  system 
employing  *  single  burst  for  detection  of  a  constant  cross  section  target  and  a  Swerling 
II  target,  respectively.  The  figures  show  probability  of  detection  as  a  function  of  the 
mean  signal-to-noise  ratio  per  burst.  It  is  useful  to  compare  these  figures  at  various 
values  of  detection  probability.  Generally  a  single-burst  probability  of  detection  of 
about  0.5  or  greater  is  needed  for  successful  tracking.  For  example,  from  Figure  7  at 
a  value  of  probability  of  detection  of  0.7,  the  loss  in  detection  sensitivity  relative  to 
undisturbed  conditions  is  about  3.5  dB  at  a  value  of  54  of  0.5  and  about  11  dB  at  a 
value  of  5<  of  1.0.  Figure  8  shews  similar  losses  in  sensitivity  relative  to  undisturbed 
conditions  for  a  Swerling  II  target  model.  With  the  exception  of  very  low  SNR, 
radar  performance  is  better  for  the  constant  target  cross  section  model  relative  to  the 
Swerling  II  model  due  to  the  decreased  variability  in  received  signal  power. 

In  Figures  9-27  the  effect  of  M  out  of  N  burst  combining  is  investigated  for 
various  values  of  M  and  N .  In  these  figures  the  abscissa  is  the  signal-to-noise  ratio 
per  look  defined  as  N  times  the  signal-to-noise  ratio  per  burst.  Thus  the  curves  are 
plotted  so  that  various  combinations  of  M  and  N  may  be  compared  on  the  basis  of 
an  equal  amount  of  transmitted  energy  per  look,  on  the  assumption  that  all  N  bursts 
are  transmitted  on  each  look. 

Figures  9-10  give  radar  detection  performance  for  the  two  target  models  for 
the  case  of  an  undisturbed  propagation  channel.  Figure  9  shows  the  effect  of  using 
multiple  bursts  per  look  against  a  constant  cross  section  target  with  no  propagation 
channel  fading.  The  optimum  number  of  bursts  per  look  is  one  and  there  is  no 
advantage  from  burst  combining.  When  the  target  statistic  obey  the  Swerling  II 
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Figure  7.  Probability  of  detecting  a  constant  target  with  a  single  burst 
for  various  scintillation  levels. 
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Figure  8.  Probability  of  detecting  a  Swerling  II  target  with  a  single  burst 
for  various  scintillation  levels. 
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model,  however,  the  received  signal  power  varies  independently  from  burst  to  burst, 
and  diversity  combining  is  effective  as  shown  in  Figure  10.  At  a  value  of  probability 
of  detection  of  0.8,  1  out  of  2  burst  combining  improves  detection  sensitivity  about  1 
dB  relative  to  no  combining  (i.e.,  1  out  of  1).  The  results  presented  in  Figures  9-10 
are  well  known.  However,  these  figures  are  included  to  allow  the  reader  to  visually 
compare  the  results  for  detection  in  a  benign  propagation  environment  with  the  results 
to  be  presented  for  a  fading  channel. 

In  Figures  11-18  quantitative  results  for  probability  of  detection  are  given 
that  show  the  effect  of  various  levels  of  scintillation  severity  on  double  threshold 
detection  for  both  target  models  of  interest  here.  For  these  results,  it  is  assumed 
that  the  signal  contribution  due  to  channel  fading  is  independent  from  burst  to  burst 
within  a  look.  This  case  of  the  propagation  channel  disturbance  is  referred  to  as  Type 
II  fading,  not  unlike  a  Sweriing  II  target  model. 

Figures  11-12  show  radar  detection  performance  for  the  two  target  models 
for  ar.  value  of  0.5  on  the  one-way  propagation  path.  This  value  of  £<  is  generally 
accepted  as  marking  the  transition  from  weak  to  strong  scattering.  Figure  11,  for  a 
constant  target,  shows  that  diversity  combining  does  offer  a  slight  advantage  for  1 
out  of  2  combining  for  detection  probabilities  greater  than  0.7.  At  lower  detection 
probability  or  with  the  use  of  more  than  two  bursts  per  look,  diversity  combining 
gives  a  loss  in  detection  sensitivity.  In  the  case  of  Sweriing  II  target  statistics,  1  out 
of  N  combining  is  seen  to  offer  some  improvement  for  all  values  of  N  shown  (2,4,8) 
in  Figure  12  when  the  detection  probability  is  greater  than  0.6.  For  example,  at  a 
probability  of  detection  of  0.75,  there  is  a  gain  of  about  2  dB  for  the  case  of  1  out  of 
2,  4,  or  8  relative  to  single-burst  detection. 

Figures  13-14  show  detection  performance  for  a  value  of  Si  of  0.707.  From 
Figure  13,  for  a  constant  target,  it  is  seen  that  there  is  a  modest  gain  available  for 
SNR  values  greater  than  about  13  dB  per  look,  with  the  amount  of  gain  depending 
on  the  number  of  bursts  per  look  and  on  the  SNR.  At  a  probability  of  detection  of 
0.75  there  is  about  1.5  to  2.0  dB  gain  in  detection  sensitivity  for  multiple  bursts  per 
look. 


From  Figure  14,  for  a  Sweriing  II  target,  a  substantial  gain  in  detection 
sensitivity  is  possible  for  values  of  probability  of  detection  greater  than  0.5.  For 
example,  at  a  probability  of  detection  of  0.75  there  is  a  4.5  dB  gain  in  using  1  out  of 
4  (or  8)  rather  than  using  a  single  burst. 

Performance  in  worst-case  scintillation,  characterized  by  an  S4  scintillation 
index  of  1.0  on  the  one-way  propagation  path,  is  considered  in  Figures  15-16.  Figure  15 
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Figure  9.  Probability  of  detecting  a  constant  target  with  no  channel 
fading  using  1  out  of  N  combining. 


Figure  10.  Probability  of  detecting  a  Swerling  II  target  with  no  channel 
fading  using  1  out  of  N  combining. 
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Figure  11.  Probability  of  detecting  a  constant  target  with  moderate  (S4 
=  0.5)  Type  II  fading  using  1  out  of  N  combining. 


Figure  12.  Probability  of  detecting  a  Swerling  II  target  with  moderate 
(Si  =  0.5)  Type  II  fading  using  1  out  of  N  combining. 
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Figure  13.  Probability  of  detecting  a  constant  target  with  Type  II  fading 
and  S4  =  0.707  using  1  out  of  N  combining. 


Figure  14.  Probability  of  detecting  a  Swerling  II  target  with  Type  II 
fading  and  S4  —  0.707  using  1  out  of  N  combining. 


Figure  15.  Probability  of  detecting  a  constant  target  with  severe  Ty  pe  II 
fading  (S<  =  1.0)  using  1  out  of  N  combining. 


Figure  10.  Probability  of  detecting  a  Swerling  II  target  with  severe  Type 
II  fading  (S<  -  1.0)  using  I  out  of  /V  combining. 
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illustrates  that,  even  for  a  constant  target,  diversity  combining  offers  a  performance 
improvement.  An  improvement  of  nearly  7  dB  in  detection  sensitivity  is  available 
when  the  receiver  operates  at  a  detection  probability  of  0.75  using  1  out  of  8  (or  4) 
detection.  For  the  case  of  the  Swerling  II  target  model  as  shown  in  Figure  16,  the 
use  of  1  out  of  8  (or  4)  combining  gives  an  improvement  of  about  4  dB  at  a  detection 
probability  of  0.5,  and  for  a  detection  probability  of  0.75  the  improvement  is  9.5  dll. 

A  comparison  of  Figures  11,13  and  15  shows  that,  for  a  constant  target  and 
Type  II  fading,  it  is  advantageous  to  use  1  out  of  N  detection,  provided  that  the  value 
of  SA  is  0.707  or  greater.  For  less  severe  scintillation,  single-burst  detection  exhibits 
better  performance.  In  the  following,  M  out  of  N  is  considered  further  for  values  of 
A/  other  than  1. 

Figures  9-16  show  the  effect  of  variation  of  the  number  of  bursts  per  look, 
A',  for  1  out  of  N  combining,  on  detection  performance  for  constant,  and  Swerling  II 
targets.  Various  levels  of  Type  II  fading  ranging  from  no  fading  to  worst  case  fading 
(S^  =  1)  arc  considered.  These  figures  illustrate  that  some  protection  against  Type 
II  fading  is  available  using  diversity  combining.  As  an  additional  example  illustrating 
the  above  point,  compare  Figures  17-18  with  Figures  7-8,  respectively.  Figure  17 
gives  probability  of  detection  using  1  out  of  8  combining  for  a  constant  target  model 
^>as  a  parametric  function  of  the  St  scintillation  index.  Results  in  Figure  7  are  for  an 
identical  case  except  only  a  single  burst  is  used.  Note  detection  performance  is  not 
strongly  affected  by  fading  when  1  out  of  8  combining  is  used.  Also  note  that  for 
a  detection  probability  of  0.7,  1  out  of  8  combining  gives  better  performance  in  a 
strong  fading  environment  than  in  a  benign  environment.  Furthermore,  performance 
is  rather  similar  throughout  the  strong  scattering  regime.  Of  course,  a  large  price 
is  paid  for  using  1  out  of  8  combining  to  detect  a  constant  target  in  an  undisturbed 
environment.  Figures  8  and  18,  for  the  case  of  a  Swerling  II  target  model,  illustrate 
a  similar  situation.  At  a  value  of  probability  of  detection  of  0.7,  a  sensitivity  loss  of  1 
dB  is  incurred  when  using  1  out  of  8  combining  as  the  propagation  environment  goes 
from  no  fading  to  worst  case  fading.  However,  the  cost  for  using  1  out  of  8  combining 
in  an  undisturbed  environment  is  less  than  2  dB  relative  to  single  burst  detection  at 
a  detection  probability  of  0.7  for  the  Swerling  II  model. 


5.1  TYPE  I  FADING. 


The  previous  discussion  has  demonstrated  the  effectiveness  of  M  out  of  N 
combining  against  Type  II  fading,  where  the  effects  of  propagation  channel  scintillation 
are  independent  from  hurst  to  burst.  If,  however,  the  time  and  frequency  scheduling 
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combining. 


Figure  18.  Probability  of  detecting  a  Swc-rling  II  target  for  various 
scintillation  levels  during  Type  II  fading  using  1  out  of  8 
combining. 
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of  bursts  is  such  that  the  effect  of  the  propagation  channel  is  identical  from  burst  to 
burst  (i.e.,  Type  I  fading),  then  diversity  combining  is  no  longer  effective.  Figures  19- 
20  show  the  effect  of  Type  I  fading  on  1  and  2  out  of  8  combining  for  values  of  S4  of  0, 
0.5,  0.707  and  1.0  for  a  constant  and  a  Swerling  II  target,  respectively.  A  comparison 
of  these  curves  to  those  of  Figures  7-8  which  give  single-burst  detection  performance 
shows  that  there  is  generally  a  loss  suffered  from  combining  in  a  Type  I  scintillation 
environment  relative  to  using  a  single  burst,  for  levels  of  scintillation  severity  where 
S4  is  greater  than  0.5. 

5.2  GENERAL  BURST  COHERENCE. 


To  investigate  the  effects  of  partial  correlation  on  the  performance  of  a 
frequency  hopping  radar,  Monte  Carlo  simulations  are  performed  using  a  Gaussian 
correlation  function  to  describe  the  second-order  statistics  of  the  received  power.  The 
radar  of  interest  is  assumed  to  change  frequency  pseudo-randomly  from  burst  to  burst. 
In  addition,  the  hopping  sequence  is  assumed  to  be  distributed  uniformly  over  some 
hopping  bandwidth,  Bhop- 

In  general,  the  received  power  decorrelates  over  time  and  frequency  as  dis¬ 
cussed  in  Section  2.4.  This  investigation  will  consider  the  case  where  the  time  separa¬ 
tion  between  bursts  is  small  enough  relative  to  the  decorrelation  time  to  assume  the 
received  signal  power  does  not  decorrelate  over  time  within  a  look.  This  assumption 
allows  use  of  a  coherence  function  that  is  a  function  of  frequency  only. 

For  this  application  a  simple  Gaussian  correlation  is  chosen  to  represent  the 
correlation  of  power  from  burst  to  burst. 


C..(/)  =  exp{-(///J)2}  (37) 

where  fi  is  the  decorrelation  frequency.  This  simple  Gaussian  correlation  function  is 
chosen  as  the  first  application.  By  generating  N  uniform  random  burst  frequencies 
over  the  hopping  bandwidth,  Bhop ,  and  using  Equation  37  to  determine  how  the  bursts 
are  correlated,  the  methods  of  Section  4.9  may  be  used  to  generate  realizations  of  the 
received  signal  power  via  Monte  Carlo  Simulation. 

Detection  performance  results  for  1  out  of  8  combining  for  a  Swerling  II 
target  are  shown  in  Figures  21-23.  The  three  figures  give  results  for  values  of  Bhopl /<* 
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Figure  19.  Probability  of  detecting  a  constant  target  for  various 
scintillation  levels  during  Type  I  fading  using  1  and  2  out 
of  8  combining. 
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Figure  20.  Probability  of  detecting  a  Swerling  II  target  for  various 
scintillation  levels  during  Type  I  fading  using  1  and  2  out 
of  8  combining. 
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of  10,  5,  and  1.  As  Bhop/ Id  increases,  the  expected  correlation  of  bursts  will  also 
decrease.  In  general,  performance  will  improve  for  larger  values  of  Bhop/ fd- 

Figure  21  gives  the  detection  performance  of  1  out  of  8  combining  for  a 
Swerling  II  target  for  various  values  of  scintillation  index  and  Bhop/ ft  =  10.  Compar¬ 
ison  of  the  curves  for  S4  =  1.0  in  Figure  21  and  Figure  13  shows  about  a  1.5  dB  loss  in 
sensitivity  at  a  probability  of  detection  of  0.7  for  Bhop/ Id  —  10  when  compared  to  the 
predicted  performance  assuming  Type  II  fading.  A  similar  comparison  of  Figure  22 
with  Figure  13  shows  about  a  3  dB  loss  in  sensitivity  when  Bhop/ Id  =  5. 

As  Bkop! Id  is  lowered  to  1,  the  preceding  comparison,  now  made  using  Fig¬ 
ure  23,  indicates  a  loss  in  sensitivity  of  about  7  dB.  More  importantly,  comparison 
of  Figure  23  with  Figure  3  (single  burst  detection)  shows  that  for  S<  =  1.0  and  a 
detection  probability  of  0.7,  1  out  of  8  combining  provides  no  diversity  gain  at  this 
raiio  of  Bhop  I  Id- 

Additional  results  (not  shown)  indicate  performance  is  essentially  identical 
to  that  expected  under  the  Type  I  fading  assumption  when  Bhop/ Id  <0.1. 


5.3  M  OTHER  THAN  1. 


Figures  24-25  show  the  effect  of  variation  in  M  on  detection  sensitivity  for 
the  case  of  an  undisturbed  propagation  channel  with  a  constant  cross  section  target 
and  a  Swerling  II  target,  respectively.  For  a  constant  target,  detection  performance 
is  best  using  4  (or  5)  out  of  8,  with  higher  values  of  M  yielding  decreased  detection 
sensitivity.  Higher  values  of  M  are  not  shown  in  the  figure  for  clarity.  For  a  Swerling 
II  target,  2  (or  3)  out  of  8  provide  the  best  detection  performance. 

Figures  26-27  show  probability  of  detection  during  Type  II  fading  of  different 
levels  of  severity  using  2  out  of  8  combining  for  a  constant  target  and  a  Swerling  II 
target,  respectively.  These  figures  can  be  directly  compared  to  Figures  17-18  which 
give  detection  performance  using  1  out  of  8  combining.  From  the  comparison  it  is 
seer,  that  1  out  of  8  yields  better  detection  performance  than  2  out  of  8  during  strong 
scintillation  conditions  (S4  >  0.5),  but  the  reverse  is  true  during  weaker  scattering 
(S4  <  0.5).  Note,  however,  this  advantage  is  dependent  on  the  assumption  above  of 
Type  II  fading,  as  may  be  seen  by  comparing  Figures  21-23  for  different  amounts  of 
burst-to-burst  correlation. 
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Figure  21. 
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Probability  of  detecting  a  Swerling  II  target  using  1  out  of  8 
for  various  scintillation  levels  assuming  Gaussian  correlation 
of  bursts,  Bkop/fd  =  10. 


Figure  22.  Probability  of  detecting  a  Swerling  II  target  using  1  out  of  8 
for  various  scintillation  levels  assuming  Gaussian  correlation 
of  bursts,  Bhop/ fd  -  5. 
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Figure  23.  Probability  of  detecting  a  Swerling  II  target  using  1  out  of  8 
for  various  scintillation  levels  assuming  Gaussian  correlation 
of  bursts,  Bh0pl fi  -  1- 
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Figure  24.  Probability  of  detecting  a  constant  target  with  no  fading 
M  out  of  8  combining. 
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Figure  25. 
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Probability  of  detecting  a  Swerling  II  target  with  no  fading 
using  M  out  of  8  combining. 


Figure  26.  Probability  of  detecting  a  constant  target  during  Type  II 
fading  using  2  out  of  8  combining. 


Figure  27.  Probability  of  detecting  a  Swerling  II  target  during  Type  II 
fading  using  2  out  of  8  combining. 
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5.1  NONCOHERENT  INTEGRATION. 


Results  are  presented  in  this  section  for  detection  using  noncoherent  integra¬ 
tion  of  bursts  for  comparison  with  the  previous  results  for  double-threshold  detection. 
Double  threshold  detection  is  sometimes  referred  to  as  binary  integration  since  it  can 
be  considered  a  modification  of  noncoherent  integration.  The  modification  is  that  the 
output  of  the  quadrature  detector  is  A/D  converted  using  1  bit  of  resolution  (thresh¬ 
old  detection)  prior  to  integration.  Since,  clearly,  some  information  is  lost  when  using 
such  a  coarse  measure  of  the  quadrature  detector  output,  one  can  expect  a  penalty 
to  be  exacted  in  detection  sensitivity  for  using  double-threshold  detection.  Results  in 
this  section  illustrate  this  loss  in  detection  sensitivity  for  Swerling  II  targets. 

Figure  28  shows  the  probability  of  detection  for  noncoherent  integration  for 
no  fading  and  various  values  of  JV,  the  number  of  bursts  per  look.  This  figure  can 
be  compared  with  Figure  10  for  double-threshold  detection  using  1  out  of  N.  From 
Figure  28  it  is  seen  that  for  noncoherent  integration  and  a  Swerling  11  target,  the 
more  bursts  per  look  the  better,  provided  that  the  signal-to-noise  ratio  is  sufficiently 
high  to  maintain  the  probability  of  detection  at  about  0.5  or  greater  and  thereby 
avoid  the  abrupt  degradation  in  detection  performance  with  decreasing  SNR.  On  the 
other  hand,  from  Figure  10,  for  double  threshold  detection  using  1  out  of  N  in  an 
undisturbed  propagation  environment,  only  1  out  of  2  is  generally  useful. 

Figures  29-33  give  results  for  the  probability  of  detecting  a  Swerling  II  target 
using  noncoherent  integration  of  8  bursts  per  look  under  various  conditions  of  scin¬ 
tillation  severity  and  degrees  of  correlation  of  the  bursts  of  a  look.  In  all  five  figures, 
curves  are  shown  for  the  case  of  no  fading  and  for  values  of  one-way  S4  of  0.5,  0.707, 
and  1.0.  Figure  29,  for  Type  II  fading,  may  be  compared  directly  with  Figures  18  and 
27  for  double  threshold  combining  using  1  or  2  out  of  8.  In  the  case  of  S4  equal  to  or 
greater  than  0.5,  1  out  of  8  gives  better  performance  than  2  out  of  8.  However,  it  is 
seen  from  a  comparison  of  Figures  18  and  29  that  noncoherent  integration  of  8  bursts 
enjoys  an  advantage  of  about  1.5  dB  in  comparison  to  double  threshold  detection 
using  8  bursts  per  look. 

Figure  30,  31,  and  32  show  detection  performance  for  the  case  of  partially 
correlated  bursts  in  a  look  with  values  of  Bhop/ fd  =  10,  5,  and  1,  respectively,  and 
may  be  directly  compared  to  Figures  21,  22,  and  23.  These  comparisons  indicate  that 
noncoherent  detection  using  8  bursts  per  look  enjoys  an  advantage  of  about  2  dB 
relative  to  double  threshold  detection  using  1  out  of  8. 
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Figure  28.  Probability  of  detecting  a  Swerling  II  target  using 
noncoherent  integration  of  bursts  with  no  fading. 


Figure  29.  Probability  of  detecting  a  Swerling  II  target  using 
noncoherent  integration  of  8  bursts  during  Type  II  Fading. 
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Figure  30.  Probability  of  detecting  a  Swerling  II  target  using 
noncoherent  integration  of  8  bursts  for  various  scintillation 
levels  assuming  Gaussian  correlation  of  bursts,  Bhop/ fd  =  10. 


Figure  31.  Probability  of  detecting  a  Swerling  II  target  using 
noncoherent  integration  of  8  bursts  for  various  scintillation 
levels  assuming  Gaussian  correlation  of  bursts,  Bhop/  fd  —  b. 
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Figure  32.  Probability  of  detecting  a  Swerling  II  target  using 
noncoherent  integration  of  8  bursts  for  various  scintillation 
levels  assuming  Gaussian  correlation  of  bursts,  Ihop.'fd  !• 
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Figure  33.  Probability  of  detecting  a  Swerling  II  target  using 
noncoherent  integration  of  8  bursts  during  Type  I  Fading. 
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Figure  33  shows  detection  performance  using  8  bursts  per  look  in  a  Type  1 
fading  environment  where  the  propagation  effects  on  all  bursts  in  a  look  are  perfectly 
correlated.  A  comparison  of  this  figure  with  Figure  20  again  gives  an  advantage  of 
about  2  dB  relative  to  double  threshold  detection  using  8  bursts  per  look. 
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SECTION  6 


CONCLUSIONS 


Severe  fading  has  a  strong  effect  on  the  target  detection  performance  of  a 
space  based  radar.  In  this  report  both  analytic  and  simulation  results  for  the  perfor¬ 
mance  of  double  threshold  detection  are  compared  to  that  of  noncoherent  integration 
of  the  returns  from  multiple  bursts  per  look.  Quantitative  graphical  results  are  given 
for  various  levels  of  scintillation  severity  ranging  from  no  scintillation  to  worst  case 
Rayleigh  fading  on  the  one-way  propagation  path.  In  addition,  the  effects  of  various 
amounts  of  correlation  between  the  bursts  in  a  look  are  considered  through  the  use 
of  Monte  Carlo  simulation  techniques.  For  this  work  a  simple  Gaussian  correlation 
function  is  used  to  specify  frequency  correlation  of  bursts  for  a  frequency  hopping 
radar.  However,  the  simulation  allows  for  an  arbitrary  correlation  in  frequency  and 
in  time  between  the  radar  returns  from  bursts  in  a  look  as  long  as  the  ‘slow  fading’ 
assumption  is  not  violated. 

For  a  nonfluctuating  target  and  worst-case  fading,  some  gain  is  available 
through  the  M  out  of  N  combining  process.  However,  as  the  level  of  scintillation 
decreases  to  less  severe  fluctuations,  burst  combining  should  be  avoided  since  better 
performance  is  possible  with  a  single  burst.  For  a  Swerling  II  target,  both  double 
threshold  detection  and  noncoherent  integration  of  multiple  bursts  per  look  offer  some 
mitigation  of  the  effects  of  fading,  depending  upon  the  severity  of  the  fluctuations  and 
on  the  correlation  properties  of  the  bursts  of  a  look. 

For  the  case  of  greatest  interest  here,  the  combining  of  8  bursts  to  form 
a  look,  there  is  a  gain  of  about  1.5  to  2  dB  in  detection  sensitivity  for  noncoherent 
integration  in  comparison  to  double  threshold  detection.  This  conclusion  is  based  on 
detection  performance  alone  and  does  not  consider  any  burst  scheduling  advantages 
of  M  out  of  N  detection. 

In  order  to  apply  the  results  for  detection  performance  given  in  this  report  to 
help  design  a  space  based  radar,  it  is  necessary  to  obtain  the  value  of  the  scintillation 
index  as  a  function  of  radar  transmission  frequency,  geometry,  ionospheric  conditions 
and  time  of  day.  The  WBMOD  computer  code  ( Fremouw  and  Secan ,  1984]  contains 
the  only  world-wide  model  that  gives  this  result.  Work  currently  in  progress  at  MRC 
involves  the  incorporation  of  a  model  of  the  probability  density  function  for  electron 
density  fluctuat;ons  into  WBMOD.  This  probability  density  function  is  necessary  to 
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answer  questions  regarding  the  percent  of  time  that  the  scintillation  index  is  above  or 
below  any  specific  value  and  is,  therefore,  necessary  to  measure  radar  performance. 
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PHYSICA1  PESEftP^H  INC 
ATTN  H  FITZ 
ATTN:  WHEUSER 

PHYSICAL  RESEARCH.  INC 
ATTN  R  DELIBERIS 
ATTN  T  STEPHENS 


TOYON  RESEARCH  CORP 
ATTN:  J  ISE 

TRW  INC 

ATTN:  H  CULVER 

TRW  SPACE  &  DEFENSE  SYSTEMS 
ATTN:  DM  LAYTON 


PHYSICAL  RESEARCH.  INC 
ATTN  J  DEVORE 
ATTN  J  THOMPSON 
ATTN:  WSCHLUETER 

PHYSICS  INTERNATIONAL  CO 
ATTN:  C  GILMAN 


USER  SYSTEMS,  INC 

ATTN  S  W  MCuANDLESS,  JR 

UTAH  STATE  UNIVERSITY 
ATTN:  K  BAKER 
ATTN.  L JENSEN 


R  &  D  ASSOCIATES 

ATTN  CGREIFINGER 
ATTN  F  GILMORE 
ATTN:  G  HOYT 
ATTN:  M  GANTSWEG 

RAND  CORP 

ATTN:  C  CRAIN 
ATTN:  EBEDROZIAN 

RAND  CORP 

ATTN  B  BENNETT 


VISIDYNE,  INC 

ATTN:  J CARPENTER 

FOREIGN 

F0A2 

ATTN:  BSJOHOLM 

FOA  3 

ATTN  TKARLSSON 

DIRECTORY  OF  OTHER 


RJO  ENTERPRISES/POET  FAC  BOSTON  UNIVERS'TY 

ATTN  A  ALEXANDER  ATTN;  MiuHAEL  MENDILLO 

ATTN.  W  BURNS 


Dist-4 


